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In diffusive Josephson junctions the phase-difference <j> between the superconductors strongly influences the 
spectroscopic features of the layer separating them. The observation of a uniform minigap and its phase mod¬ 
ulation were only recently experimentally reported, demonstrating agreement with theoretical predictions up to 
now - a vanishing minigap at <j> = n. Remarkably, we find that in the presence of intrinsic spin-orbit coupling 
a giant proximity effect due to spin-triplet Cooper pairs can develop at (j> = it, in complete contrast to the 
suppressed proximity effect without spin-orbit coupling. We here report a combined numerical and analytical 
study of this effect, proving its presence solely based on symmetry arguments, which makes it independent of 
the specific parameters used in experiments. We show that the spectroscopic signature of the triplets is present 
throughout the entire ferromagnetic layer. Our finding offers a new way to artificially create, control and isolate 
spin-triplet superconductivity. 


I. INTRODUCTION 

Spin-polarized superconductivity in which the Cooper pairs 
reside in a so-called triplet state is currently attracting much 
attention, topical examples including chiral p-wave pairing, 
singlet/triplet mixing in superconductors without an inversion 
centePS intrinsic coexistence of ferromagnetism and super¬ 
conductivity in heavy-fermion compounds®! and lately topo¬ 
logical superconductivity^ In all of these cases, triplet 
Cooper pairs emerge which offer the tantalizing prospect of 
combining superconductivity and spintronics'P 

Spinless singlet Cooper pairs can penetrate into other ma¬ 
terials via tunneling from a superconductor across their mu¬ 
tual interface, a phenomenon know as the proximity effect 11 . 
When the adjacent material is a ferromagnet, the penetration 
distance depends on the texture of the magnetic exchange 
field, which may convert the singlet pairs into triplets. If 
triplet components with spin projection parallel to the field are 
generated, these correlations can penetrate much further into 
the ferromagnet, being known as the long-range triplet (LRT) 
component 12 . The primary route to creating the LRT com¬ 
ponent has so far been to use magnetically inhomogeneous 
structures 12 13 or non-equilibrium distribution functions and 
intrinsic triplet superconductor^. However, it was recently 
established that spin-orbit (SO) coupling can also provide the 
necessary rotatiorPEHto create LRT. Experimentally, such SO 
coupling can be provided either by the structural asymmetry 
in thin-film hybrid structures 17 , giving rise to generic interfa¬ 
cial spin-currentJ-L] or due to the intrinsic crystal structure, 
i.e. provided by noncentrosymmetric materials^. 

The superconducting proximity effect is a phase-coherent 
phenomenon that can be probed in e.g. Josephson junctions 
where the density of states is highly sensitive^ to the su¬ 
perconducting phase difference (ft. Experimentally, Ref. El 
reported measurements for Josephson junctions with a nor¬ 
mal metal (SNS), that were consistent with the theoretical 
prediction 22 that the density of states evolves from a finite 
minigap due to the superconducting correlations (</> = 0) to 
the absence of such a minigap (</> = 7r). This can be under¬ 
stood intuitively: the proximity effect would be suppressed 
when the order parameter in each superconductor is equal in 


magnitude but opposite in sign, resulting in superconducting 
correlations “averaging” to zero and a featureless density of 
states in the center of the normal-metal region. 

We here discover that when SO coupling is present in a 
magnetic Josephson junction, the opposite effect takes place. 
Remarkably, the SO coupling in the junction instead gives rise 
to a giant, triplet-induced proximity effect at (ft = n, as shown 
schematically in Fig. [T] We prove that the reason for this is 
that the SO coupling forces the triplet Cooper pairs to have 
the opposite parity symmetry compared to the singlet Cooper 
pairs with respect to the center of the junction. When <fi = i r, 
the singlet correlations are antisymmetric across the junction 
whereas the triplets are symmetric, meaning that the proxim¬ 
ity effect survives even in the centre of the junction and is 
solely due to LRT Cooper pairs. In other words the phase dif¬ 
ference, which is an experimentally tunable quantity, can be 
used to remove the spin-singlets and keep only triplets even 
with a homogeneous exchange field. Previous attempts to sep¬ 
arate spin-polarized Cooper pairs from the singlet component 
have required magnetic inhomogeneities, which are experi¬ 
mentally challenging to control, so the inclusion of SO cou¬ 
pling for this purpose represents a significant step forwards. 
We prove analytically, solely from symmetry arguments, that 
the giant triplet proximity effect always occurs near the centre 
of the junction, making its existence independent of the spe¬ 
cific junction parameters used. Whereas the density of states 
can be probed locally, as can spin injection and tunneling in 
superconductor^, we also show numerically that the gener¬ 
ation of spin-polarized Cooper pairs for 7r-biased Josephson 
junctions in fact persists throughout the entire system. The 
fact that our prediction is based solely on symmetry and that 
its spectroscopic fingerprint is independent of where the local 
density of states is measured in the system makes it a very 
robust effect. 

The remainder of the article will be organised as follows: In 
Section [II] we introduce the quasiclassical theory and notation 
for the trilayer Josephson junctions and discuss the analytic 
consequences of the inclusion of SO coupling in this formal¬ 
ism. In Section[III]we present numerical results demonstrating 
the emergence of the giant triplet proximity effect at </> = 7r, 
and moreover show that this persists throughout the interstitial 
junction layer. Finally, we conclude in Section[IV]with a brief 
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(a) Standard SNS junction 

DoS vs. E 



(b) SFS with spin-orbit coupling 

DoS vs. E 



FIG. 1: The proximity effect in the standard SNS Josephson junction 
(a) manifests in a minigap that closes with increasing phase differ¬ 
ence (f> between the superconductors, closing entirely for <j> = n. 
Conversely, the proximity effect in the SFS junction with intrinsic 
SO coupling in the junction direction (b) results in a giant peak in 
the density of states at zero energy for <j> = n. The different be¬ 
haviour is due to the symmetry of the anomalous Green’s function: 
for SNS. the singlet component is symmetric around the centre of the 
junction at f = 0 and antisymmetric for <j> = n whereas the SFS case 
retains a symmetric triplet component at f = n. For the examples 
used in-text we take the junction direction to lie in the ^-direction 
and use bulk values for the superconductors. 


summary and discussion. 


II. THEORY 


It is instructive to briefly consider the behavior of the prox¬ 
imity effect in an SFS junction without SO coupling as a func¬ 
tion of the phase difference. In this case the quasiclassical 
Usadel equation 24 in the linearized regime without SO cou¬ 
pling reads D F d^f± + 2ie±f± = 0, where e± = e ± h z 
and f± = ft ± f s for energy e, magnetisation exchange field 
h, diffusion constant Dp in the ferromagnet and singlet and 
triplet anomalous Green’s functions f s and f t respectively. 
The Usadel equation describes the diffusion of the condensate 
into the adjacent material, and the corresponding Kupriyanov- 
Lukichev boundary conditions 25 to the superconducting in¬ 
terfaces take the form £ Lpd z f± = T/Bcse 1 ^ at z = 0 
and C L F d z f ± = ±/ BC se ,0R at z = L F where / BC s is the 
bulk, Bardeen-Cooper-Schrieffer anomalous Green’s function 


in the superconductors, L F is the length of the ferromagnet 
and Q is the interface parameter. <p F and Or denote the re¬ 
spective superconducting phases. The solution in the middle 
of the junction reads: 


f± 


±/bcs cos(k±L/2) , i 4 , n , 

{Lism(k±L) [ + 


( 1 ) 


where k± = \J2ie± j D F is the wavenumber. By direct in¬ 
sertion, one observes that when the phase difference & = 
4>l — ( Pr = 7r, the superconducting proximity effect vanishes 
completely since J± = 0. This holds for all energies and re¬ 
gardless of whether h = 0 or h 0. Since this takes place 
at <j> = 7r both in SNS and SFS junctions, one might be led 
to think that this is a robust phenomenon. However, we now 
show that in the presence of spin-orbit interactions, this is no 
longer the case. Not only is there a strong superconducting 
proximity effect present in the junction at = tt, but for cer¬ 
tain cases it is in fact the maximum triplet proximity effect that 
can be obtained in the junction. This is in stark contrast to the 
conventional picture of a vanishing minigap at <f = 7r reported 
previously for SNS junctions-: 

The giant triplet proximity effect can be established solely 
on symmetry arguments, making it independent of the specific 
junction parameters employed in an experiment. To see this, 
we have derived the Usadel equations in terms of a Ricatti- 
parametrization 23 making it very suitable for numerical cal¬ 
culations, which we use later in this paper. The Usadel equa¬ 
tion with SO coupling expressed in terms of the full Green’s 
function was derived in Ref. [T33 — in its Ricatti-parametrized 
form in the ferromagnet, one obtains (see Ref. [27] for our full 
derivation): 


D F [pi') + 2(dk'y)N')(dk'y) S j 
= — 2 *e 7 — th.(a 7 — 70 ;*) 


+D F 


AA 7 — 7 H* A* 4- 2 (Ay + 7 H* )N (A* + 7 A 7 ) 


+2 iDp 


{d k '))N(A* k +jA k ')) + (A k +'yA* k j)N(d k 'y) 


■( 2 ) 


Here A = {A x , A y , A z ) is the SO field with components 
Ak = A^a a , using the summation convention over re¬ 
peated indices with a = x,y, z, and index k indicates an 
arbitrary choice of direction in Cartesian coordinates. The 
vectors h and a are the corresponding three-component ex¬ 
change field and vector of Pauli matrices. The 7 , 7 are the 
Ricatti-parametrised matrices characterising the quasiclassi¬ 
cal Green’s function g R : 


zR = (N(I+ 77) 2AT 7 \ 

9 V —2./V7 - 7 V(/ + 7 7 y’ 


(3) 


with normalisation matrices N = (I — 77 ) _1 and N = 
(/ — 77 ) ” 1 and identity matrix /. The 7 operation denotes 
complex conjugation and e —) (—e). Similarly, the boundary 
conditions become: 

<9fc7i = -r~r{I- 7i 72 )^ 2(72 - 7i) + + *7r4lfc, 

dkl2 = y— ~{I - 727r)-^i(72 - 7i) + *^fc72 + *72(4) 

^2S,2 
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In the last line of Eq. 0 we see that the first order deriva¬ 
tive couples to the component of the SO field in the corre¬ 
sponding direction. This is the key requirement for the ap¬ 
pearance of the giant triplet proximity effect at <f> = n. Con¬ 
sider therefore an SFS junction oriented along the z-direction 
as in Figure |T] and a choice of SO coupling vector aligned 
perpendicular to the interfaces ( A z ^ 0). For concrete¬ 
ness and to give more transparent analytical results, we set 
A = (0,0 ,aa x — aa y ). This choice corresponds e.g. to 
Rashba-type coupling with broken inversion symmetry in the 
n = x + y direction, with the transverse motion of the elec¬ 
trons restricted, effectively corresponding to motion in a one¬ 
dimensional wire. We underline that the non-vanishing triplet 
proximity effect in 7r-biased Josephson junctions predicted 
here survives even if one were to include components of the 
SO coupling field {A x , A y } / 0 and if one were to only in¬ 
clude the a x or <i y term in A z above. We can then write down 
the linearised Usadel equations in the limit of weak proximity 
effect, where we have \^\ij <C 1, N ss 1 in the ferromagnet 
such that 


2 7 = ( fa f* + ft] (5 ) 

7 \ft~fs hi )■ { ’ 

The origin and main features of the giant triplet proximity 
effect can now be identified analytically by considering the 
low-energy regime e = 0 and setting the exchange field to 
h = h z z, where the equations become: 

{d 2 z - 4 a 2 )f a <j + 4<ra(l - ai)d z f t - Uaa 2 f- a - a = 0, 
D F d 2 z f s + 2 ih z f t = 0, 

D F d 2 z f t + 2 ih z f s - 8 D F a 2 f t 

+ 2D F a(l - %)d z f u - 2D F a(l + i)d z f n = 0. (6) 

with a =t) i- When there is zero phase difference between 
the superconductors, the anomalous Green’s function for the 
singlet pairs, f s , is a symmetric function with respect to the 
middle of the junction. This can be seen immediately from 
the general form of the solution of f s and the boundary condi¬ 
tions, and is equivalent to what happens for conventional SNS 
and SFS junctions. When the phase difference is equal to n, 
however, f s (and thus necessarily its second derivative) is an¬ 
tisymmetric. Now, it follows from Eq. [2] that when the SO 
coupling has a component in the junction direction it neces¬ 
sarily introduces a first-order derivative term. Performing the 
operation z —> (—z) on Eq. 0, this means that the function 
subject to the first order differential must be a symmetric func¬ 
tion, provided it is not constant. In the linearised regime 0 
we can see explicitly what this entails: the functions and 
hi must be symmetric around the middle of the junction, and 
it is clear that a nonzero component of the anomalous Green’s 
function will remain at zero energy even in the 7r-biased junc¬ 
tion. Since the density of states in the ferromagnet is given 
by 


D(e) = TV {N(I + 77)}A (7) 


such that the expression at zero energy in the linearised regime 


becomes 


25(0) = 1- 


l / s ( 0)| 2 

2 


\.m \ 2 

2 


l/n(Q)l 2 

4 


l/u(0)l 2 

4 


,( 8 ) 


an experimental signature of this effect would be an enhanced 
zero-energy density of states due to the triplets. Without SO 
coupling, this typically occurs for exchange fields fulfilling 
the resonant conditioners] h ~ E g , where E g is the minigap 
energy. However, the mechanism for the enhanced density of 
states due to triplets in the present work is fundamentally dif¬ 
ferent from previous literature as it originates from a parity 
effect due to SO coupling which is present independently on 
the exact junction parameters. We underline that the purpose 
of Eq. ([8]) is to illustrate that the triplets will enhance the den¬ 
sity of stated 30 ! whereas in the actual numerical calculations 
shown in the figures we have solved the non-linearized Us¬ 
adel equation and computed the density of states via Eqn. 0 - 
Although the precise phase difference that produces maximal 
peak at zero energy depends on the strengths of the exchange 
field and SO coupling (and we leave a more detailed exposi¬ 
tion of the model dependencies for later work) the <j> = p case 
remains significant and can also be maximal. We emphasize 
again that this discovery is the complete opposite of what has 
been thought to be the case up to now in Josephson junctions, 
namely a vanishing superconducting proximity effect when¬ 
ever the phase difference is n. 


III. RESULTS 


In what follows, we present for the first time a solution to 
the Usadel equation with spin-orbit coupling in the full prox¬ 
imity regime. In Fig. [2] we show an example of the spectro¬ 
scopic profile for varying phase difference between the su¬ 
perconductors in an SFS junction with SO coupling, high¬ 
lighting the generation of a zero-energy peak in the density 
of states at (/) = 7T. We choose an in-plane exchange field 
h = 10Ay for ease of experimental application, a bulk super¬ 
conducting coherence length = 30 nm and SO coupling 
strength a = 0A/Lp, i.e. normalised to the ferromagnet 
length Lp, here chosen to be 15 nm such that the relevant 
quantity Lp/^s = 0.5. Qualitatively similar behaviour is ob¬ 
served for most other choices of exchange field orientation 
and SO coupling strength. A comparison with the standard 
SNS and SFS junctions without SO coupling is provided and 
the giant proximity effect at <j> = n is immediately clear. Since 
the singlet component f s vanishes when (/) = ", the remain¬ 
ing features are entirely due to the triplets and in this case en¬ 
tirely due to the LRT component. In effect, the result reported 
here serves as a way to fully isolate the triplet correlations 
regardless of the junction parameters in 7r-biased Josephson 
junctions. We note that quantitatively, even when h > A, 
the proximity effect and resulting enhancement of the density 
of states here is much larger than experiments measuring the 
same quantity for superconductor-ferromagnet hybrids with¬ 
out SO coupling 31 3 21, where the deviation from the normal- 
state is about 1%. 








4 



cp = 0 - <p = 0.25jt - <p = 0.5jr .cp = 0.7571 - (p = n 


FIG. 2: The local density of states D(e) at z = Lf /2 for phase 
difference <j> between the superconductors of an SFS junction with 
spin-orbit coupling aligned in the junction direction, with exchange 
field h = 10Ay and spin-orbit coupling a = OA/Lp . The giant 
triplet proximity effect at tj> = 7r manifests as a large peak in the D(e) 
at e = 0. A comparison with the standard SNS and SFS junctions 
without SO coupling is also provided. 


Another important question pertains to how the manifesta¬ 
tion of the giant triplet proximity effect depends on the dis¬ 
tance from the superconductors and also on asymmetries be¬ 
tween the two interfaces due to e.g. different transparencies. 
Remarkably, the effect is virtually independent of the distance 
from the superconducting interfaces: the spectroscopic peak 
originating from the presence of spin-polarized Cooper pairs 
persists all the way up to the interfaces and hardly changes 
throughout the junction as shown in Fig. [3] For reference, the 
spatial distribution of the anomalous Green’s function is pro¬ 
vided in Fig. [6] in the Appendix. For substantially stronger 
spin-orbit coupling, a spatially oscillatory contribution from 
the short range triplet component with spin projection perpen¬ 
dicular to the field can also arise as one moves toward the 
interfaces, but the proximity effect still remains completely 
dominated by triplets throughout the entire system, i.e. not 
only near the middle of the junction (see Figs.[5]an More¬ 
over, we have also checked numerically (not shown) that the 
spatial dependence remains unchanged even for asymmetric 
junctions where one interface is as much as twice as transpar¬ 
ent as the other, i.e. a ratio of barrier parameters C 1 /C 2 = 2. 
This suggests that the predicted effect should be very robust 
and facilitates its experimental observation. 

In Figure[4]we show the intricate dependence of the density 
of states of the 7r-biased junction on the strengths of exchange 
field h = hAy and SO coupling aLp, which is highly non¬ 
monotonic in a. We see that as the field strength increases a 
more narrow spectrum of SO coupling will generate a giant 
peak at zero energy, with the optimal SO coupling decreasing 
slightly for higher field strengths. 

In Figure[2]we saw that the peak at zero energy was present 
for all phases, increasing significantly as </>—>■ n. This is the 
predominant behaviour observed for the majority of the a- 
h parameter space, with a few exceptions. In particular, we 
note that when the exchange field is oriented in the junction 
direction, i.e. in the same direction as the SO coupling, in- 
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FIG. 3: (Color online) The spatial dependence of the local density 
of states D(e, z) for a phase differences <j> = n between the super¬ 
conductors of an SFS junction with spin-orbit coupling. Parameter 
values are the same as in Fig. [2] As seen, the giant proximity ef¬ 
fect is virtually independent on the position from the interfaces, thus 
persisting throughout the entire system. Inset: top view of the DoS 
surface plot. 



FIG. 4: The local density of states D(e = 0) at z = Lf/2 for a 
7r-biased SFS junction with spin-orbit coupling vector aligned in the 
junction direction (corresponding to a nanowire setup) as a function 
of magnetization exchange field h = hAy and strength of spin-orbit 
coupling olLf aligned in the junction direction. 


creasing the SO coupling can take the density of states profile 
from having a zero-energy peak for all phase differences as in 
Fig. 0 to a profile with a qualitative change from minigap at 
0 = 0 to peak at <f> = n, as shown in Fig. [5] The reason for the 
appearance of a full minigap in Fig. [5] at zero phase difference 
in spite of the presence of an exchange field can be understood 
by noting that for large spin-orbit coupling, the triplet compo¬ 
nent becomes suppressed whereas the singlet one remains un¬ 
affected. In this case the SO coupling stabilizes the minigap, 
which now persists throughout the interstitial junction layer, 
and leads to a magnetically tunable minigap 27 . However, the 
isolation of spin-polarized LRT Cooper pairs at = tt sur¬ 
vives, providing an enhanced density of states D(e = 0) > 1, 
indicating that it is a robust phenomenon. 
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FIG. 5: The local density of states D(e) at z = Lf/2 for phase 
differences <f> between the superconductors of an SFS junction with 
spin-orbit coupling aligned in the junction direction, with exchange 
field h = 3 A z and spin-orbit coupling a = 2/Lf. At small phase 
differences D(t) displays a minigap while retaining a zero-energy 
peak at </> = 7r. 

IV. DISCUSSION AND SUMMARY 

In the present analysis we restricted the form of the SO cou¬ 
pling vector to lie along the junction direction, as in the case 
of one-dimensional wires, although we have also verified that 
the giant proximity effect is maintained for weak transverse 
SO coupling in the xy-plane. A common alternative way of 
introducing intrinsic SO coupling compared with using a non- 
centrosymmetric crystal is to use the interfacial asymmetry of 
thin-film junctions, and in this case Rashba-Dresselhaus cou¬ 
pling will lie in the try-plane. This case will be explored in 
more detail in Ref. [27] We also considered here the coexis¬ 
tence of the exchange field and the SO coupling, which may 
also be taken as a model of the case where these features ex¬ 
ist in separate, adjoining layers, such that the SO coupling 
is induced, e.g. by deposition of a heavy transition metal or 
compound. Such a setup could be easier to integrate into cur¬ 
rent devices, and the consequences for supercurrent genera¬ 
tion will be explored in the future^. 

In summary, we have shown that the inclusion of intrin¬ 
sic spin-orbit coupling in the ferromagnet of an SFS junction 
can result in a pure spin-triplet state without any singlet su¬ 
perconductivity, and that this state can persist throughout the 
ferromagnet. Remarkably, this occurs for a phase difference 
of 7r between the superconductors, dispelling the commonly 
held view that the proximity effect should be suppressed at 
(f> = tt. In fact, in many cases the peak at zero energy will 
be significantly enhanced compared with the features at zero 
phase difference, displaying what we term the giant proximity 


effect in 7r-biased junctions. The inclusion of intrinsic spin- 
orbit coupling is therefore a very promising resource for har¬ 
nessing triplet superconductivity for improved functionality in 
spintronic devices. 
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Appendix A: Spatial variation of anomalous Green’s function 

To examine the spatial variation of the singlet and triplet 
components of the anomalous Green’s function /, we write 

f = ifs + d ■ o)io y . (Al) 

Now the singlet component f s transforms as a scalar under 
spin rotations while the triplet component transforms as a 
vectoi^EHand is described by the d-vector d = ( d x , d y , d z ) = 

((/u - /tt)> ~*(/w + /tt)> 1 2 3 4 /t)/ 2 - In Fi g- |6] we show how 
the magnitude of the real part of the singlet and triplet compo¬ 
nents vary throughout the ferromagnet for the examples dis¬ 
cussed in the main text (the imaginary component is zero at 
e = 0). With an in-plane exchange field h = lOAy as shown 
in Fig. |6la), the LRT component is given by the dominant 
Re (d x ) [throughout the entire length of the ferromagnet. The 
short-ranged triplet component is defined as the component of 
the d-vector parallel to the exchange held, and it is evident 
that in this case the contribution from the short-ranged triplet 
increases towards the superconductor interfaces but remains 
minimal compared with the LRT component. 

Increasing the strength of SO coupling can give rise to os¬ 
cillations between the long- and short-ranged components, as 
is evident in Fig. |6j&), with an out-of-plane exchange held 
h = 8A£. In this case both d x and d y are symmetric and show 
the long-ranged component with a maximal peak around the 
middle of the junction. The short-ranged component d z is zero 
in the centre of the junction and increases correspondingly as 
the LRT decays. In both cases the singlet contribution is neg¬ 
ligible, such that the triplet component dominates throughout 
the junction. The strength of the exchange held has a minimal 
effect on the frequency of oscillation in comparison with the 
strength of SO coupling. 
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